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1. Introduction 

Ultrarelativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the 
Large Hadron Collider (LHC) are able to reach temperatures high enough to create and study 
the quark-gluon plasma (QGP) in a controlled experimental environment. So far, fluid-dynamical 
descriptions that take into account only the effects of shear viscosity have been considered sufficient 
to describe the time evolution of the QGP created in such collisions. 

However, QCD is a nonconformal field theory and, in principle, one should not neglect the 
effects of bulk viscosity on the dynamics of the QGP. The bulk viscosity is expected to be maximum 
near the (pseudo) QCD phase transition and, if it is sufficiently large, it can drive the effective 
pressure of the system to become negative [1, 2], driving fluid dynamics out of its domain of 
validity and changing the freeze-out process of the system. This effect was referred to as cavitation 
in Ref. [2]. 

In this contribution we discuss the onset of cavitation in event-by-event hydrodynamic simula¬ 
tions of heavy ion collisions at LHC energies. We find that if the bulk viscosity coefficient becomes 
large even in a small temperature domain, 170 <T < 190 MeV, it can already considerably affect 
the evolution of the system. We further estimate the value of bulk viscosity, around this narrow 
temperature region, that is required to make the effective pressure become negative. 

2. Model 

The time evolution of the hot QCD matter produced is solved using relativistic dissipative fluid 
dynamics. The main equations of motion are the continuity equation for the energy-momentum 
tensor, =0, where = eu^u'^ — A^''{Po+U) + 7i^'', withebeing the energy density, 

Pq the thermodynamic pressure, the fluid velocity, H the bulk viscous pressure, and the 
shear-stress tensor. Here, we introduced the projection operator A^'' = onto the 3-space 

orthogonal to the fluid velocity. The equation of state, Po(£)^ reflects chemical equilibrium and is 
taken from Ref. [3], corresponding to a parametrization of a lattice QCD calculation matched to a 
hadron resonance gas calculation at lower temperatures. We assume that the baryon number density 
and diffusion are zero at all space-time points and our metric convention is = diag(+l, — 1 — 
1 - 1 ). 

The time-evolution equations satisfied by H and are derived from kinetic theory [5, 6] and 
solved numerically using the MUSIC hydrodynamic simulation [7, 8]. We solve 

Tnfl + n = —^0 — ^nnllP-|-An;c7r^'"o'^v , (2.1) 

pn^'' =2rio^'' -b A^rnlla^''. (2.2) 

For the sake of simplicity, most of our transport coefficients are calculated using the Boltzmann 
equation near the conformal limit [6]. For the shear viscosity, we use the temperature dependent 
coefficient (HQ-HH parametrization) defined in Ref. [9]. For the bulk viscosity, we consider two 
parametrizations: the bulk viscosity defined in Refs. [10, 11], which corresponds to a parametriza¬ 
tion of calculations from [12], for the QGP phase, and [13] for the hadronic phase. These two calcu¬ 
lations are matched around Tc = 180 MeV and the value of at this temperature is (^/s{Tc) ~ 0.3. 
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The other choice of bulk viscosity is a parametrization of calculations from lattice QCD [2], for the 
QGP phase, and, once more, [13] for the hadronic phase. In this case, the calculations are matched 
around Tc = 180 MeV with a much larger value of ^/s, ^/s{Tc) « 1. So both parametrizations of 
i^/s are the same in the hadronic phase and are very similar in the QGP phase (at high temperatures, 
both choices of 1^/5 are basically zero). However, they differ considerably in the phase transition 
region, for temperatures 170 < T <190 MeV. In the following we will see that, even though the 
viscosities are only different in such narrow temperature domains, they are still able to considerably 
modify the time-evolution of the system. Both parametrizations are shown in Fig. 1(a). 



Figure 1: (Color online) The left panel shows the bulk viscosity over entropy density ratio as a function of 
temperature (normalized by Tc = 180 MeV). The solid blue line corresponds to parametrization 1 and the 
dashed red line to parametrization 2. The right panel shows the initial temperature profile employed in all 
simulations discussed in this work. 


As a test case, we consider the hydrodynamic simulation of one perfectly central (impact pa¬ 
rameter equal to zero) Pb-Pb collision event at LHC energies. The initial state considered is calcu¬ 
lated using the MC-Glauber model, with the initial energy density of the system being proportional 
to the density of participants. The initial time is taken to be Tq = 0-6 fm and we further assume that 
the system starts at rest and in local thermodynamic equilibrium, i.e., u‘ (Tq) = IT (tq) = (fo) = 
0. The initial entropy of the medium is determined in such a way that the multiplicity of charged 
hadrons measured at the LHC is described. The initial temperature profile is plotted in Fig. 1(b). 


3. Results 

In Figs. 2(a) and (b) we show the space-time profile of H/Pq in the chosen hydrodynamic 
event for the two parametrizations of i^/s shown in Fig. 1. The H/Pq profiles are shown in the 
r = plane, along the x = y axis. 

We see that the effective pressure of the system, Peff = Pq -|- H, can become very low in some 
spacetime points for both parametrizations of i^/s considered. Such small values of pressure occur 
exactly in the phase transition region, where the bulk viscosity has a peak. The second parametriza¬ 
tion, which has the largest bulk viscosity around the phase transition region, leads to the formation 
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Figure 2: (Color online) Spacetime evolution of H/Fo in a hydrodynamic simulation of a central Pb-Pb 
collision. The left panel shows the simulation that employed parametrization 1 of ^/s while the right panel 
shows the same simulation using parametrization 2. 


1=2.6 fm 


1=4.6 fm 


1=6.6 fm 



Figure 3: (Color online) Spatial profiles of temperature (upper panels) and H/Pq (lower panels) for three 
time steps of our simulation, T = 2.6 fm (left panels), T = 4.6 fm (middle panels), and T = 6.6 fm (right 
panels). 


4 













Onset of cavitation in ultrarelativistic heavy ion collisions 


Gabriel S. Denicol 


of a negative pressure hypersurfaee that eneloses the whole system. That is, it is impossible for 
a fluid element to reaeh the eonstant temperature freezeout hypersurfaee (for all freezeout tem¬ 
peratures, Tfo <180 MeV) without passing through a negative pressure region. However, for 
parametrization 1, where (^js < 0.3, even though the effeetive pressure ean reaeh small values, 
it is still positive in most spaee-time regions. We have eheeked that, by inereasing the value of bulk 
viseosity around Tc in parametrization 1 from ^ /s{Tc) = 0.3 to ^ /s{Tc) = 0.6, eavitation starts to 
oeeur. 

In Figs. 3 we show the transverse spatial profile of the temperature and bulk viseous pressure 
over thermodynamie pressure for three different time steps, z = 2.6 fm (left panels), z = 4.6 fm 
(middle panels), and z = 6.6 fm (right panels). The upper panels show the temperature profiles 
while the lower panels display the 11/Pq profiles. This time evolution was ealeulated with the 
seeond parametrization of l//s, whieh has the largest value around l//s{Tc) ~ 1. We already 
showed in Fig. 2 that regions with negative effeetive pressure will oeeur in the fluid when sueh 
large values of bulk viseosity are employed. 

We see that, at the initial times, there is always a ring of negative effeetive pressure formed 
near the edges of the system. These rings of negative pressure gradually propagate towards the 
eenter of the fluid and, at T = 6.6, the eenter of the fireball is already filled by a negative pressure 
domain. In the temperature profiles shown in the upper panels of Fig. 3, we see the formation of 
"holes" or poekets of small temperatures ereated in the system due to the negative effeetive pressure 
aehieved in these regions. Also, we see that the bulk viseosity is so large that it eonsiderably reduees 
the expansion rate of the system, leading to a very uniform temperature profile in the eenter of the 
fireball. These effeets resembles what happens in a first order phase transition, where the expansion 
of the system eonsiderably slows down in the mixed phase region. 

4. Conclusions 

In this work, we investigated the onset of eavitation in event-by-event hydrodynamie simu¬ 
lations of heavy ion eollisions. We showed that, if the bulk viseosity is of the order of (^/s 1 

around the QCD phase transition region, the fluid-dynamieal evolution of the system will generate 
negative effeetive pressures in eertain regions of the fluid. This behavior indieates the breakdown 
of the fluid-dynamieal deseription and, in the eontext of hydrodynamie simulations of heavy ion 
eollisions, will lead to ehanges in the way the freezeout proeess of the system is developed. We 
have also found that, if the bulk viseosity is ss 0.3 around the QCD phase transition region, 
eavitation will not happen, but the bulk viseosity will have signifieant effeets on the evolution of 
the system. 

This work was supported by the Natural Seienees and Engineering Researeh Couneil of Canada. 
G. S. Denieol aeknowledges support through a Banting Fellowship of the Natural Seienees and En¬ 
gineering Researeh Couneil of Canada. 
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